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ABSTRACT 


This paper studies axially symmetric vibrations of a liquid-filled 
poroelastic thin cylinder saturated with two immiscible liquids of 
infinite extent that is surrounded by an inviscid elastic liquid. By 
considering the stress free boundaries, the frequency equation is 
obtained. Particular case, namely, liquid-filled poroelastic cylinder 
saturated with single liquid is discussed. When the wavenumber is 
large, the frequency equation is reduced to that of Rayleigh-type 
surface wave at the plane boundary of a poroelastic half-space. In 
this case, the asymptotic expressions of Bessel functions and 
modified Bessel functions are used. In both general and particular 
cases, the case of the propagation of Rayleigh waves in a poroelastic 
half-space is obtained. The parameter values of Columbia fine 
sandy loam saturated with air-water mixture are used for the 
numerical evaluation. In all the cases, phase velocity as a function 
of wavenumber is computed and presented graphically. From the 
numerical results, some inferences are drawn. 
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1 INTRODUCTION 


N general, wave propagation problems in liquid-saturated porous media is attracting many researchers due to its 

importance in various diversified areas of Science and Engineering. These problems have been studied 
extensively in Soil Mechanics, Seismology, Acoustics, Structural Engineering, Geophysics, and Biomechanics. Pipe 
lines in the oil exploration and water supply systems may be surrounded by some kind of liquids. These pipe lines 
are cylindrical and filled with liquids. Study of wave propagation in liquid-filled cylindrical pipes surrounded by 
liquid provides some information about the strength of pipe lines since wave characteristics depend on the strength 
of the material. Biot developed constitutive equations and the equations of motion in a poroelastic solid saturated 
with a single fluid [1]. This well established theory could not address the issues in important domains, namely, 
Petroleum Engineering and Bone Mechanics, wherein porous solid is saturated with two or more fluids. To this 
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extent, Tuncay and Corapcioglu developed a theory of wave propagation in a poroelastic solid saturated by two 
immiscible Newtonian fluids in the Eulerian framework by employing the volume averaging technique [2]. A study 
of body waves in poroelastic solid saturated by two immiscible Newtonian fluids is presented [3]. In the paper [3], 
the existence of three compressional waves and one rotational wave in an infinite porous medium is shown 
analytically. Santos et al. proposed a method to determine the elastic constants for isotropic porous media saturated 
by two fluids [4]. Sahay et al. developed a set of equations to describe the low frequency seismic phenomenon in 
porous solid which is inhomogeneous and anisotropic by using the volume average technique [5]. Hunyga 
developed a general dynamical model for porous media saturated by two immiscible fluids [6]. Lo et al. derived a 
general set of coupled partial differential equations to describe dilatational wave propagation through an poroelastic 
medium saturated by two immiscible fluids [7]. Propagation of axisymmetric waves through compressible, inviscid 
fluid contained in a cylindrical, elastic shell is studied by Lin [8]. Using exact three dimensional equations of linear 
elasticity, the frequency equation for the vibrations of a fluid-filled circular cylindrical shell is derived by Ram 
Kumar[9]. Sharma and Gogna solved the problem of elastic wave propagation in a cylindrical bore in poroelastic 
solid and derived the frequency equation for empty and fluid-filled bores [10]. Axisymmetric wave propagation 
along fluid-loaded cylindrical shells is investigated both theoretically and experimentally by Plona et al. [11]. The 
longitudinal waves in homogeneous anisotropic cylindrical bars immersed in a fluid are studied by Dayal [12]. The 
laws of propagation of axisymmetric normal modes in a hollow cylinder filled with and surrounded by fluid media 
are investigated [13]. Vashishth and Khurana studied the propagation of wave along a cylindrical bore hole 
embedded in an anisotropic porous solid saturated by viscous fluid [14]. Harmonic radiation from a spherical surface 
vibrating at the center of a fluid-filled circular cylindrical cavity is investigated [15]. The propagation of elastic 
waves along a cylindrical bore hole filled with and without fluid embedded in an infinite porous medium saturated 
by two immiscible fluids is studied by Arora and Tomar [16]. In the paper [16], the dispersion equation of Rayleigh- 
type surface waves along the boundary of a poroelastic half space saturated by two immiscible fluids is obtained. 
Rayleigh waves in poroelastic half space and a porous solid lying over an elastic solid are examined by Tajuddin et. 
al. [17,18]. In both the papers, it is assumed that porous solid saturated with a single fluid. A comparative dispersive 
study between the bone without marrow and bone with marrow is made in the framework of Biot's theory of 
isotropic poroelasticity by Malla Reddy et al. [19]. Axially symmetric vibrations of fluid-filled and empty 
poroelastic circular cylindrical shell of infinite extent are investigated by Ahmed Shah [20]. Flexural vibrations of 
poroelastic cylinder of various thicknesses immersed in a fluid are studied in the paper [21]. Vibrations in a fluid 
loaded poroelastic hollow cylinder surrounded by a fluid under plane strain conditions are investigated by Shanker 
et al. [22]. wave propagation in composite hollow sphere and fluid-filled sphere surrounded by fluid are analyzed in 
the papers [23, 24]. Radial vibrations of poroelastic cylindrical shell immersed in an acoustic medium is studied in 
the paper [25]. 

In the present paper, axially symmetric vibrations of a liquid-filled poroelastic thin cylinder of infinite extent 
saturated with two immiscible liquids that is surrounded by an inviscid elastic fluid are investigated. Propagation of 
waves in a liquid-filled poroelastic thin cylinder saturated with a single liquid is recovered as a particular case. In 
both the cases, the frequency equation of Rayleigh wave is obtained. The rest of the paper is organized as follows. In 
section 2, governing equations and solution of the problem are presented. Boundary conditions and frequency 
equation are given in section 3. Particular cases are presented in section 3.1. Numerical results are discussed in 
section 4. Finally, conclusion is given in section 5. 


2 GOVERNING EQUATIONS AND SOLUTION OF THE PROBLEM 


Let (r,0,z)be the cylindrical polar coordinates. Consider a homogeneous, isotropic, infinite poroelastic thin 
cylinder saturated with two immiscible liquids loaded internally and externally by inviscid liquids. The axis of the 
cylinder is in the direction of z- axis as shown in Fig. 1. The equations of motion in the absence of body forces for 
low frequency wave propagation in a poroelastic solid saturated with two immiscible fluids are given under [16]: 


1 > > > > > > > > Ou 
V((a, +p WV -u,+a,V-u,+a,;,V-u,)V-(G,Vu, )+c (viy, ) +e (VV) = P, aie y 
> > > > > oO a 8 = (1) 
> > > > > u 
V(a V u, +4, V -U,+a,,V 0) -c (V-V) = A, a? » V(a V u, +ayV-U,+a,,V-U,)—C,(V,—V,) = Py age > 
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4,=Ky, a, =a, =K,a,S,(A,+K,)/D, a; =a, =K,a,(1—-S, (A, +K,)/D, 


A 
an =K S7-a, XK, EA a» =A, =K KS (A-5, )(l—-a@, )/D, 


1 


A 
a,, =K,(1-S,} (1-a, XK, Tee iD, D=K,(1-S,)+4,+K,S, 
1 


ci =(1-a,) S? 4 /KK,,, ©, =(l-@,)°d-S,) 4, /KK,,, S,=a,/(1-a@,), 
d * * 
A, =PrS,(1-S;,), Pras A —P>), 


u. and u, are displacement vectors in the poroelastic solid and in the liquid of phase i respectively. v. and v, are 
velocity of solid phase and in the liquid of phase i respectively; p, is the volume averaged density of porous solid, 
p;, L; and K, are the volume averaged density, viscosity, bulk modulus of fluid phase i respectively. œ, is the 
volume fraction of the solid, œ, is the volume fraction for the liquid phase i respectively; K is the intrinsic 
permeability of solid and K,, is the relative permeability of liquid phase i. G,, is the shear modulus of solid and 


K,,, is the frame or drained bulk modulus, Pp, and p} are the pressures of non-wetting liquid and wetting liquid. 


Fig.1 
Liquid-filled thin cylinder surrounded by a liquid. 


The stress in solid and the pressure in liquids are: 
> > > > > T 2 > 
(c,)=(@,V-u,+4,V-u,+a,V-u,)/ +G, (Vu, +(Vu, y — zV% D), 
(z) =(a,,V-u,+a,,V-u,t+a,,V-u, I, (2) 


(t,) =(a,,V-u,+a,,V-u,+a,,V-u, I, 


> > 
where J is unit tensor matrix. For axially symmetric vibrations, let u, =(u,,0,w,), u,=(,,0,w,), and 


u, =(u,,0,w,) be the displacement vectors of solid, first liquid and second fluid, respectively; which are functions 


of r,z and t. The pertinent components are obtained from the field equations Eq. (1) and are given under: 
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u, =—(Ap,K\ (pr) +4,p.K (pr) + AsikK (psr )e E”, 

w, =(4,ikK (pyr) +A ikK (p,r) +4 CEK (pst) pK prye, 

u, =-(4,8 p K (pyr) +4,8 pK (pyr) + A563 pK (pr) +4,8p K (pyr) -mm4 ikK (pr E”, 

w, (AS PK (D7) + AD DK P) AED Kp AS PK Aat)-maM ALEK (pr) K pe, 6) 
u, =-(4;P;K (pr) +4,p,K (pyr) -momp A;ikK (pre E”, 

w, = (APK (pyr) +A pK (pat) ZA EK P) -pK (pre E”, 


where @ is the frequency of the wave, k is the wavenumber, t is time, i is the complex unity, 4,,4,,A,,A4, and 4; 


are all constants, K „(x ) is the modified Bessel functions of second kind of order n, and 


1 1 
P, =k l-4), j =1,2,3,4,5, éi =¢7((ayE, +43) + (0E, +43)" —40'm,,dE,)*)/ 2ayE,, 


1 1 
E =° (E, +45) - (0E, +04) -40m0 E,)?)/ 20E, 65 =° (gE y+ (gE) -40M paE,))/ Udy, 


1 2 2 
> m m 
é =c° (aE, —((ayE 4) —40'm,dyE,)?)/2a,E,, E =c*(m,, a= =), 
My My» 
4 ee 4) 
E =p) (a; toO /@ (maz =Ma»), E, =M TAM / mnan =M 24»), 


= 2 — 
E, =al; -anan / O (Mya MAn), E, =AyM3/ (Mya — M.A), 
2 2 
ô, =E% +E,, l =1,2, ô, =E,ć +E,, l =3,4, 
i (ci +c,) ic, ic, ic, ic, 


my s > My = > My = > My = Py > M3 = Py — $ 
(72) (7) (72) (72) (72) 


In Eq. (4), c is the phase velocity. Here poroelastic thin cylinder under consideration is loaded and surrounded by 
two different liquids. Therefore, we shall proceed to find the stresses of a liquid. For axially symmetric vibrations, 
the displacement potential function ¢ of a liquid satisfies the following equation: 


Ob 109 Fb _ 1 OY 


ôr? rôr @z? V? at? (5) 


The displacement components of inner liquid u, = (u, 0,w ,) are obtained from the solution of Eq. (5) and the 


equations Uy = $ and Wy = x, 
r z 


The displacement components are 


uy =Ba,L@ re!” wy =B ik qyr eo, (6) 


zu 

where B, is constant, /,(x) is the modified Bessel function of first kind of order n, and q, = kd. The 
if 

inner liquid pressure is 


2 


Py =P pg = pB I (qyr E. (7) 
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> 


The displacement components of outer liquid u, = (u, 0,w p ) are obtained on similar lines, which are given 


oj 


under: 


Uy = Body AY Gy rye we) Wop = B,ikH ® ayr a (8) 


21 
where B, is constant, H® (x) is the Hankel function of first kind of order n, and q, y =k a-2. The outer 
of 
liquid pressure is 


Po = Por a B,J, oy rie? ou (9) 


By substituting the displacement components from Eq. (3) into Eq. (2), the solid stresses and liquid pressures are 
obtained, which are 


T, +T, +T, =(4,M,,(r)+4,M,, (7) +4,M (7) +4,M ar) +A M (rye, 
7, =(4M a(r) +4M y(r) +4M >, ()+4A,M yr) +4; Mp (rJe EPL 
T, =(4,M,,(r)+4,M,,(r)+4,M,,(r)+4,M y rei E, 
T, =(4M „ (r) +A,M ,(r)+4,M p) +A, M prp E, 


(10) 


where, 


M (r) =G, +a; a, Ha, + + (ap +yy +4y))5; p? RK On AGREN, 
i r 
M (r) = (a; +4, +4 +(a +ay 4 Ay Os (Ps -k°))K (pr), M (r) =2G IKK (psr), 
M (r)=-2G,;pikK (pr) M(r)=M,(r)=0, M„(r)= (k? —p;)K (psr), 
M, (r) = (a, +a, Xp? -k DK, (pir), M,, (r)= (4,3 +4,,0; (Ds -k DK, (pr), 
M a(r) = (a; +a, ô; p? =k DK, (pr), M (7) = (a; +a, 8 p3 =k DK (pr), 


M (r), M (r) are similar expressions as M, (r), M, (r) with p,,6; and p,,6; are replaced by p,,6; and 
P4,0,, respectively. 

M „(r) is similar expressions as M „(r) with p,,6, are replaced by p,,ô;, respectively. 

M ,.(r), M,,(r) are similar expressions as M, (r), M,,(r) with pp ô and p,,0; are replaced by p,,ô, and 
P4,0,, respectively. 
2 
2 


M p (r), M,,(r) are similar expressions as M (r), M4,(r) with p,,0, and p,,6, are replaced by p,,6, and 


P4,0,, respectively. 


3 BOUNDARY CONDITIONS AND FREQUENCY EQUATION 


For perfect contact of solid and liquid at the interface, the displacements and stresses must be continuous. Thus the 
boundary conditions to be satisfied on the curved surface r =a to be stress free are turned to be as follows: 


Tr +7, ET, =—p 


IZ 


y> Tmt +T, =—Py>s T, =9, 


q =0, T, =0, u,—u, =0, u, —u y =0. (11) 
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Substitution of Eqs. (6)-(10) into Eq. (11) result in a system of seven homogeneous equations in seven unknowns 
A,,A,,A,,A,,4;,B, and B,. In order to obtain a non-trivial solution of this system, the coefficient matrix must be 
singular. Accordingly, the frequency equation is obtained, which is given by 


|4, |=0, i,j =1,2,3,...7, (12) 


y 
where, 


A, =M, j =12,345, As =p 0l qra), 4g =9, A, =M, j =1,2,3,4,5, Ax =0, Ay =p H q0), 
Ay =M, j =1,23,4,5, 4y=0, dy =0, Ay, =M, j =1,2,3,4, A =A% =4y=0, 

Ay, =M; j =1,2,3,4, Ass =4s=4=0, Aa=-PK (pa), Aa=-P:K (pa), As =e =9, 

Ass =ikK (p,a), As =-G 1G), Ag =0, 4, =4Ag, An =4Ao, Ay =Ay =9, 

Ait) =Ags, An =F Hi Gy )- 


3.1 Particular cases 


In this subsection, two particular cases when ka is large are considered. The first case is solid saturated with two 
liquids and the other case is solid saturated with a single liquid. Both the cases result in the frequency equations of 
Rayleigh waves in a poroelastic half space. 


3.1.1 Half space saturated with two liquids 


When the wave length is very small or the wavenumber is large (i. e., when ka — œ); the frequency equation Eq. 
(12) is reduced to that of Rayleigh-type surface wave at the plane boundary of a poroelastic half-space saturated 
with two immiscible liquids. In this case, the asymptotic expressions of Bessel functions and modified Bessel 
functions [26] are employed and the pertinent frequency equation Eq. (12) is obtained, which is given under: 


|=0, i,j =1,2,3....7, (13) 


where, 


2 2) 22 2) g2 
B= 2G, +(a, Brie, +a, +a, + (Ay +a + sy )O, Sy, Bis (7) = (a; +a, +53 + (a + yy +Ayy )O; G3 > 


B5 =2iG,&, B= prs By= Ube =A j =1,2,3,4,5, Bea Rape 
By =-21G,€, By =-2iG,6), By =B; =0, B; SiGe, By =B; =0, 
By =(4, +n 6, Es By =@; +a; E55 B, =By =By =0, 

B, =@ +a, )& 5 By, =(@y +a, és, B; =B; =B; =9, 


1 1 Di 
t 1 . ood 
By =(1-é'), By =(1-6)?, Bo By 0, Bs l, Bo a 72) Bar =0, 
if 


1 2 ail 

S eee 

By, =B,, l-6), j =1,2,3,4,5, B =0, B,, ATA 
of 


B,,, B,,are similar expressions as B,, ,B,, with €,,0, and &,,0, are replaced by £,,5; and &,,6;, respectively, 
B, B,, are similar expressions as B,,, B,, with é,,0, and &,,0, are replaced by €,,6, and &,,6;, respectively, 


B,,, Bare similar expressions as B,,, B,, with €,,6, andé, ô; are replaced by €,,5; and &,,6;, respectively. 
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3.1.2 Half space saturated with a single liquid 


When the presence of one liquid, out of the two immiscible liquids and outer fluid are neglected, so that 
a, =a, =a, =0, Py =V =0. Then the frequency equation Eq. (12) reduces to 


=0, i,j S123 4, (14) 


where, 


2 2 
Ca =G, +a -30s +46; (pi -KYK (pr) + G,K, (pr), 


Ci 54s Cua =Aye, Ca =4;> Cy 54n; Cy =Aps, Cy =0, Cy 545 
Cy =Ay, Cy =Cy =0, Cy =Ag, Cy 54o, Cy =0, Cu =Aoe- 


C,, are similar expressions as C, with p,,6, and p,,0; are replaced by p,,6, and p,,6;, respectively. The 


frequency equation Eq. (14) is analogous to the Eq. (3.8) of the paper [19]. when ka >, the Eq. (14) reduces to 
the frequency equation of Rayleigh wave in a fluid-filled poroelastic half space saturated with single liquid, which is 
given under 


|D,|=0 i,j =1,2,3,4, 
Dı =By, Dy =By, Da =B; Dy =B Dy =B; Dy =By, Dy =B, Dy =0, 


(15) 
D, =B D, =B D,, =B,,=0, D, =B Dp =B Dg =0, Dy =B 


51? 52? 61? 62? 


4 NUMERICAL RESULTS 


For numerical work, we consider poroelastic thin cylinder saturated with non-viscous liquids. In this case, 
L; =0, which implies c, =0,j =1,2. Frequency equations are investigated for a cylinder made up of Columbia fine 


sandy loam saturated with air-water mixture [16]. The physical parameters for the above solid are given in Table 1. 
Oil is mostly available in sandstone reservoirs. Cylindrical casing pipes containing perforations are used to pump the 
oil to the surface. These casing pipes sometimes made of sandstone itself. The velocity of sound in water (v,,) and 


kerosene (~ „) are taken to be 1.432x10°m /s and 1.3251x10°m/s , respectively. Density of water (p,) and kerosene 


(Py ) are taken to be 1000kg /m° and 820.1kg /m’, respectively. Employing these values in frequency equations, 
we obtain a implicit relation between phase velocity and wavenumber. 


Table 1 
Parameter values. 
Material parameter value 

K, 1.44x10°N / m? 
K, 35x10°N /m° 
K, 0.145x10°N / m? 
K, 2.25x10°N / m? 
Gy, 1.02x10°N / m? 
Ps 2650kg / m° 
Pı 1.1kg /m a 
Pz 997kg / m° 
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Bisection method is used to solve the frequency equations and numerical process is performed in MATLAB. 
Numerical results are depicted graphically in Figs. 2-3. Fig. 2 depicts variation of phase velocity against 
wavenumber in the case of liquid -filled poroelastic thin cylinder saturated with two liquids surrounded by a liquid. 
From Fig. 2, it is seen that the phase velocity is periodic in nature. Fig. 3 depicts variation of phase velocity against 
wavenumber in the case of liquid-filled poroelastic thin cylinder saturated with a single liquid. From Fig. 3, it is 
observed that the phase velocity is periodic in nature. when ka — œ , liquid-filled thin cylinder immersed in a liquid 
reduces to half-space separating two liquids. The phase velocity of the Rayleigh mode at the half space is computed 
and are depicted in the figures. From figures, it is seen that the phase velocity of Rayleigh wave in the case of half 
space saturated with two liquids is more than that of half space saturated with a single liquid. 


—s~— Thin cylinder 
—<— Half space 


Fig.2 

Variation of phase velocity against wavenumber in the case 
of liquid-filled poroelastic solid saturated with two liquids 
o we Sia pes a surrounded by a liquid. 


Phase velocity 


Wavenumber 


6 —— Thin cylinder 
—— Half space 


Phase velocity 


i Fig.3 

1 Variation of phase velocity against wavenumber in the case 
of liquid-filled poroelastic solid saturated with a single 
’ Eo aE IER AD pP. g Oa TAS A liquid. 


Wavenumber 


5 CONCLUSIONS 


Axially symmetric vibrations in a liquid-filled poroelastic thin cylinder saturated with two immiscible liquids 
surrounded by a liquid are investigated in the framework of Tuncay and Corapcioglu. Frequency equation in the 
case of poroelastic cylinder saturated with a single liquid is obtained as a particular case. In both the cases, Rayleigh 
waves in the poroelastic half space are investigated. Phase velocity against wavenumber is computed and depicted 
graphically. From numerical results, it is observed that the phase velocity of Rayleigh mode in the case of half space 
saturated with two liquids is more than that of half space saturated with a single liquid. 
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